Nonlinear Spectral Characterization of Discrete Data 
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The explicit analytical expression of the Nonlinear Fourier Transform (NFT) of a finite set of data 
is provided. Then a simple recursion relation for the NFT is constructed as a function of the spectral 
parameter. These tools provide a complete characterization of the nonlinear coherent structures 
(solitons, breathers, ...) present in numerical or experimental data representing the solution, at a 
given value of time, of a nonlinear evolution equation (e.g. of the nonlinear Schrodinger family). 
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Introduction. The inverse scattering (or spectral) 
transform (1ST) is a method for constructing and solving 
nonlinear evolution equations (NEE) Q. One of those, 
the nonlinear Schrodinger equation (NLS) ||, has be- 
come a paradigm as being the model for a wide vari- 
ety of physical situations. Actually the relevant equation 
generally turns out to be nonintegrable, e.g. NLS with 
forcing and damping, which renders necessary a numer- 
ical analysis ||. Then, a fundamental question is the 
characterization of the solution, as a standard Fourier 
analysis will not individualize the intrinsically nonlinear 
structures (as solitons). The same question is even more 
essential when dealing with experimental data. 

Another fundamental aspect is the question of the dis- 
cretization of an integrable continuous model, which rises 
many interesting questions such as numerical induced 
chaos [Q, instabilities in integrable systems || or round- 
off error growth ||. 

Moreover, localized excitations (solitons) in nonlinear 
chains has been a subject of intense research for their 
obvious physical interest Q || ||. It has been proved 
for instance that the discreteness generically helps energy 



localization |1C | 



LP 



wave collapse 

tence of moving localized modes 



that it is a mechanism controlling 
and also that it allows for the exis- 
ts] or standing ones 



In all these studies, the existence and properties of 
such localized excitations result merely from a direct ob- 
servation of the numerical solution of the model. In short 
a soliton is recognized from its shape, its velocity from 
the velocity of the point of maximum amplitude, the fre- 
quency of a breather is obtained from its Fourier spec- 
trum, etc... 

However there exists a tool, the nonlinear Fourier 
transform (NFT), which allows to characterize com- 
pletely the nonlinear modes. It is based on 1ST and has 
been developed numerically either for the sine-Gordon 
family |D| to modelize fluxon dynamics in Josephson 
junctions, or for the Korteveg-de Vries model [|l6| to mod- 
elize water waves. In that case, the method has been suc- 
cessfully applied to experimental data jl7| ], for instance 
to prove that solitons generically occur in wave trains. 



The drawback however is on the one side the complexity 
of the method, and on the other side the lack of analytical 
formula for the NFT. 

This letter is intended to provide two results, essential 
for the nonlinear analysis of discrete data. The first one 
is an analytical explicit formula (|ll| ) for the NFT of a dis- 
crete set of data. More precisely, taking advantage of the 
finiteness of numerical or experimental data, we obtain 
the explicit expression of the Jost solutions, and hence of 
the complete spectrum characterization, corresponding 
to a discrete finite potential. The second one is a quite 
simple recursion relation ( |l7| ) for a very straightforward 
numerical construction of the spectrum. 

Everything is done in the context of the Ablowitz- 
Ladik spectral problem jl8| (discrete version of the 
Zakharov-Shabat spectral problem), which means that 
the nonlinear Fourier analysis applies here for instance 
to the NLS family, that is 



id t q(n) = q(n+l) + q(n—l) - 2q(n) ± 

\q(n)\ 2 [q(n+l) + q{n-l)] + o(n) 



(1) 



where a(n) is any perturbation term. The method works 
without change for all nonintegrable deformations of the 
others integrable NEE related to the Ablowitz-Ladik 
spectral problem, like mKdV, Toda, self-dual network, 
KdV, sine-Gordon p9|, or discrete Stimulated Raman 
Scattering ^0|. Last, the tool can be successfully ap- 
plied to experimental data, as soon as the measurements 
are assumed to represent a train of envelope wave pulses 
which can be considered to vanish outside a finite inter- 
val. Indeed in such a case the model equation for the 
envelope is likely to obey a member of the NLS family. 

Statement of the problem. Given a finite discrete set 
of 2(L + 1) complex valued data q(n), r(n) having signif- 
icant values on n = 0, . . . , L and vanishing outside, the 
problem is to determine its nonlinear Fourier spectrum 
(the scalar case will enter as a reduction, e.g. r = ^fq 
to recover (0)). This is done by using these data as the 
potential of the Ablowitz-Ladik scattering (or spectral) 
problem: 

/in(n+l) - //n (n) = g(n+l)/i 2 i(n+l), 
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/i 2 i(n+l) - fc _1 /i2i("-) = r(n+l)/in(n+l), 
^i 2 (n+l) - fc/xi 2 (n) = g(n+l)/i 2 2(n+l), 
M22(«+l) - = r(n+l)/ii 2 (n+l), 



(2) 



for n = — 1, . . . , L — 1. A dependence on the spectral pa- 
rameter k is understood everywhere for the matrix valued 
function /z(fc,n) with elements nij. 

The Jost solutions ^ are then defined by the following 
discrete integral equations 
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q(i)n 22 (i) 



(3b) 
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Nonlinear Fourier spectrum. The spectral transform 
of {g, r} is then completely determined by the Jost solu- 
tions /i + and as follows. An essential function is the 
reflection coefficient a ± (fc) defined as 
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Due to the finiteness of the support of the potential 
{q, r}, the Jost eigenfunction, and hence a ± (fc), can be 
defined in the whole complex fc-plane. Then the knowl- 
edge of the reflection coefficient in the complex plane al- 
lows to determine completely the nonlinear Fourier spec- 
trum, which consists of a continuous part (radiation) re- 
lated to the values of on the unit circle \k\ = 1, and of 
a discrete part (solitons) constituted by the poles kj of 
or in |fc| > 1 and kj of a + in |fc| < 1. It is our purpose 
here to give an explicit formula for a in terms of the 
complex variable k. 

Explicit solution. Introducing, for (k, n) the new 
functions 
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(5) 



the integral equations (||) become explicitly solvable in 
terms of the new matrix v(k,n). This solution can be 
written in (k,n) 



'21 
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The matrix X(n) reads 

L 

X(n)= J] (1-Ai), A, 

i=n+l 

with X(L) — 1, or else 

1 '0 Q 
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(8) 



= E U ?J + E ^ <°> 

odd i=l v 7 oven i=2 v 7 

with the following definitions of the functions of the two 
variables (I, n) : 



k h q(ji)k j2 r(j 2 )...k n q(ji), 



Rl = J2 k ]1 r( n )k-^q(j 2 ) . . . k j 'r(ji) 

Ui=J2 k- jl q(ji)k h r(j 2 ) • • • fc*r(j,), 

^ = E k jl r(h)k- j2 q(j 2 ) . . . k-*qVi). (10) 

The sums run on all possible different / indices ji , ordered 
for each value of (I, n) as n + 1 < ji < j 2 < . ■ ■ < ji < L. 

Then we may compute from X(n) the reflection coef- 
ficients a ± (fc) which from their definition (0) read 
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Note: the reduction r(n) = ±q(n) leading e.g. to the 
NEE (|), has the counterpart @ a+(l/k) = ±cF(jfe). 
These two relations are indeed compatible with (p"l|). 

Data reconstruction. The data {r, g} can be recon- 
structed from their spectral transform {a + , a~} because 
the solutions ^ ± (fc,n) of (||) also solve 
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where the integral on ( runs on the unit circle. 
The system (||) at order zero in k then gives 

q(n+l) = -/J^ 2 (n) ( -~ 1 \ r(n+l) = -/4l(") 



(13) 



where fJ-i2i n Y is the coefficient of k~ x in the Laurent 
expansion for k — > oo of Hi 2 (k, n), and ^21 ( rl )^ the coef- 
ficient of k in the Taylor expansion for k — > of /i^i (fc, n). 

These formulae are useful to prove the consistency of 
the method, that is that the reconstructed data are the 
same as the starting ones (in particular that they stay 
in the class of finite support). We do not develop this 
lengthy proof here but rather make the following remark: 
the above formulae allows one to filter the data just by 
setting to zero any part of the spectrum. For instance, 
if one needs information on the underlying nonlinear co- 
herent structures, one may set to zero the contribution 



a (C)j ICI = 1 to the radiation, and then solve (12) 



Of course, doing this the reconstructed data do not stay 
in the class of finite support potentials, and this proce- 
dure will be meaningful for potentials going to zero fast 
enough at both ends. 

Recursion relation. From its very definition (^|), the 
matrix X(n) can actually be considered as a function 
also of the variable L (remember that L is the dimension 
of the data support). We shall then denote it by X L (n) 
and for instance X m (n) has to be understood as obtained 
from the data of {q(n),r(n)} for n = 0, 1, . . . , to, and zero 
for n (f. [0, m] . Then our purpose here is to derive a recur- 
sion relation for the matrix elements of X m (n), in order 
to obtain a recursion relation for the corresponding spec- 
tral data a^(fc), spectral transform of {q(n),r(ri)} for 
n = 0, 1, . . . , to. 

The spectral transform ( p"l| ) is first conveniently rewrit- 
ten as 



Pn 
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(14) 



with the following definitions 

m m 

S(m) = ]T *(n)X3(n), P(m) = 5>(n)*S(n), 

71=0 1=0 

p(k, n) = k~ n q(n), s(k, n) — k n r(n). (15) 

Then only the diagonal elements of X m (n) have to be 
considered. 

By a careful rewriting of the element X^ +1 (n) out of 
@, and by using that obviously X[\(n) = A^\ +1 (n) = 1, 
we obtain 

■m — 1 

X™+\n) = X™(n) + s(m + 1) £ p(j + l)X^(n). 



(16) 



A similar recursion relation is immediately deduced for 
S(m) defined above, and consequently for a~ . The same 
computation being made for XJ£ +1 (n) leads finally to 
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If the data to be analyzed are given on [0, L], the spectral 
transform a^(k) is constructed from the above recursion 
relation by running m from to L. To illustrate the 
preceding results we consider two examples. 

Example: two-point data. In the case of two-point 
information (L = 1), by direct application of the above 
formula, we get 



r(0) + fcr(l) 
1 + fc" 1 g(l)r(0) 



1 + Ag(0)r(l) 



(18) 



An obvious and natural consequence is that two point 
data cannot represent a pure soliton solution which would 
require that a vanish on the unit circle \k\ = 1. The 
important information is that one can readily conclude 
about the presence of a soliton if |g(l)r(0)| > 1 (which 
ensures a pole of a~ in |fc| > 1) or/and if |g(0)r(l)| > 1 
(which ensures a pole of a + in \k\ < 1). 

2nd Example: truncated soliton. We consider the one 
soliton solution (reduction r = —q) centered on m = uq 
and cut to the left of to = and to the right of to = L 
defined by (770, #0 arbitrary constants and b > 0) 



q{m) 



o irt a m+e 



cosh (bm — bno) 
sinh b 



(19) 



for m = 0, 1, • • • , L and q(m) = for m < and m > L. 
By using the recursion formula ([l7]) the spectral trans- 
form can be explicitly computed giving (£ = ke~ lria ) 
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(20) 



A simple analysis for L n shows that has a pole 
for |fc| < 1 if and only if no > 0. This means that the 
presence of a soliton in a wave train is detected with this 
method by analyzing only a finite portion of it. In the 
example considered, it is in fact enough to catch half of 
the soliton. 

Conclusion. The method is quite simple to use and it 
is of wide application: the formulae ( fll] ) provide an ex- 
plicit function of the complex variable k which allows to 
completely characterize the nonlinear Fourier spectrum 
of the data {r(n),q(n)}. In particular the roots of the 
denominator of aT in the region |fc| > 1, and those of 
a + in |fc| < 1, furnish the nonlinear coherent structures 
present in the data. For practical applications, the re- 
cursion relations ( |l7| ) provide a quite efficient numerical 
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code to generate the nonlinear Fourier spectrum of any 
finite data. These relations possess a continuous coun- 
terpart which will be studied, together with details and 
applications, in a forthcoming paper. 
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